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Abstract; 

The  ionic  rare  gas  clusters  Ar+  and  Xe+  and  rare  gas-hydrogen  clusters 

n  n 

(Ar^H)  ,  (Ar^^Hj)  ,  (Xe^H)  and  (XenH2)+  are  studied  by  the  semiempirical 

diatomics -in- ionic -systems  (DIIS)  method.  The  Ar+  clusters  (n  >  3)  are  seen 

n 

to  have  a  structure  of  a  linear  Ar+  core  surrounded  by  n-3  neutral  or  almost 

neutral  Ar  atoms.  For  Xe*  (n  >  3),  a  symmetrical  Xe^  ionic  core  with  the 

geometry  of  regular  pyramid  is  formed.  The  rare  gas-hydrogen  clusters  with 

one  H  atom  have  a  simple  R^(RH)+  structure  with  k  neutral  rare  gas  atoms 

attracted  to  the  (RH)+  molecule  by  polarization  forces.  Two  H  atoms  can  bind 

with  Ar  atoms  to  form  quasistable  clusters  Ar  which  dissociate  to 

n  / 

(n-l)Ar  +  H  +  (ArH)+  through  a  high  barrier  of  roughly  0.75  eV.  Two  H  atoms 
and  one  Xe  ion  are  shown  to  form  a  collinear  valence -bound  (XeHH)+  cluster 


□  □ 


The  simplest  small  clusters  are  formed  by  (1)  rare  gas  atoms, ^  (ii)  rare 

2-4 

gas  atoms  and  nonpolar  molecules  and  (ill)  rare  gas  atoms  and  polar 
molecules.'’  '  In  all  these  clusters,  electrons  are  strongly  localized  on  the 
component  molecules  and  rare  gas  atoms  such  that  the  particles  are  bound 
together  only  by  dispersion  (van  der  Waals)  and  in  case  (iii),  weak 
polarization  forces.  Such  clusters  are  usually  called  van  der  tfaals  (VDW) 
clusters.  The  rare  gas  atoms  can  also  form,  however,  clusters  in  which 
electrons  are  partly  delocalized,  contributing  to  some  valence  forces  between 
the  particles  involved.  In  particular,  such  delocalization  takes  place  in 
neutral  rare  gas -halogen  and  ionic  rare  gas  clusters. 

As  an  example  of  neutral  rare  gas -halogen  clusters,  we  can  mention 
XenCi.  The  ground -state  Xe^Ci  clusters  usually  seem  to  be  of  the  weakly-bound 

a 

VDW- type.  However,  afe.  initio  calculations  of  the  XeCi  complex  show  that  the 

large  electron  affinity  of  the  Ci  atom  leads  to  a  small  but  noticable  electron 

9 

delocalization  which  contributes  about  50%  of  the  XeCi  dissociation  energy. 

This  delocalization  affects  the  structure  of  the  Xe  Ci  clusters  and  is 

n 

responsible  for  the  coupling  between  their  ground  state  and  excited  ionic 

+  -  9 

(i.e.,  charge  transfer,  XenCi  )  states.  In  these  excited  states  the 

strongly-bound  quasistable  molecules  Xe^Ci  are  formed  which  have  been  studied 

10-12  13  + 

experimentally  both  in  the  gaseous  and  solid  phases.  The  Xe^Ci 

molecule  can  attract  more  Xe  atoms,  mainly  within  Xe  solid,  forming  Xe*Ci  (n 

^  „  9,14-15 

>2)  clusters. 

The  ionic  rare  gas  clusters  R+  have  been  studied  extensively  both 

experimentally  and  theoretically.  The  smallest  clusters  R*  have  been 

16-21  22  23 

detected,  in  particular,  in  Ar,  Xe  and  Kr.  initio  calculations  of 

Ne^  and  Ar^,^  ^  as  well  as  semiempirical  calculations  of  Ar^  ^  ^  and 
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+2729  + 

Xe^,  *  show  that  the  triatomic  clusters  have  a  symmetrical  linear 

structure  with  the  central  atom  bearing  roughly  50%  of  the  charge.  Clusters 

+  30 

with  more  than  three  atoms  are  found  for  all  rare  gas  atoms,  i.e..  He  , 

n 

Ne+,3*  Ar+, 27,32-36  Kr+,37  and  Xe+.38,3^  The  large  rare  gas  clusters  R+  with 
n  n  n  n  °  n 

n  >  10  demonstrate  some  magic  numbers  in  their  stabilization  energy  dependence 

on  n  33,35,37  43  «phe  structure  of  the  R+  clusters  with  n  >  3  is  not 

n 

35  + 

completely  clear.  According  to  an  experimental  study,  the  Ar^  ion  is  the 

core  in  the  Ar+  clusters.  Theoretical  Investigations  which  take  into  account 

+  26 

the  electron  delocalization  support  this  finding,  at  least  for  Arn  and 

+  29  +  + 

Xe  ,  and  show  that  R,  ions  can  also  form  a  core  of  R  (n  >  4)  clusters.  The 
n  4  n 

involvement  of  Rt  and  rT  ions  in  the  R+  cluster  formation  allows  one  to 
3  4  n 

33  41  42 

conclude  that  models  with  the  charge  localized  on  one  or  two  ’  centers 
are  unsatisfactory,  at  least  for  small  R*  clusters.  As  shown  in  Ref.  27  and 
29,  the  charge  can  be  shared  even  by  more  than  four  atoms. 

Rare  gas  atoms  can  form  not  only  the  homonuclear  clusters  R*  considered 
above  but  also  ionic  heteronuclear  clusters  which  consist  of  either  different 
rare  gas  atoms  or  rare  gas  atoms  and  valence -active  atoms.  The  charge 
distribution  and  the  structure  of  heteroatomic  ionic  clusters  depend  strongly 
on  the  relation  betwet*.  the  ionization  potentials  of  the  atoms  involved.  For 
example,  the  (XenNem)+  cluster  is  obviously  formed  by  a  Xe*  ionic  core  and 
neutral  Ne  atoms,  since  the  ionization  potential  of  He  is  much  higher  than 
that  of  Xe .  In  the  case  of  valence -active  atoms  such  as  hydrogen,  the 
situation  becomes  more  complicated  because  of  the  formation  of  the  valence - 
bonded  molecules  H^,  Hj,  and  (RH)*.  The  (Rl^)*  systems  have  been  studied 
in  the  context  of  the  R+  +  and  R  +  chemical  reactions  (see,  for  example, 
Refs.  44-46).  The  (ArHj)*  potential  energy  surface  has  been  calculated  by  the 


47 


The  study  of  the  rare  gas -hydrogen  ionic  clusters  has 


DIM  method. 


4 


concentrated  mainly  on  ArH^  clusters,^®  whose  structure,  by  analogy  with 
+  49  + 

HeH^,  has  been  suggested  to  be  a  triangle  with  Ar  at  the  vertex.  This 

structure  is  supported  by  £&  initio  calculations.^  The  ionization  potential 

of  Ar  (15.76  eV)  is  much  higher  than  that  of  H  (13.6  eV)  ,  so  that  the  charge 

is  located  on  the  H  atoms.  A  similar  charge  distribution  is  expected  for  the 

(Hj^H^  ionic  clusters,  which  have  been  studied  experimentally,^  since  the 

ionization  potential  of  (15.43  eV)  is  also  much  higher  than  the  H 

ionization  potential,  initio  calculations  of  rare  gas -hydrogen  clusters 

+  52 

with  one  H  atom  have  been  performed  for  N®nH  clusters.  Due  to  the  large 

difference  between  the  Ne  and  H  ionization  potentials,  the  charge  is  strongly 

localized  on  H+,  and  the  Ne  atoms  are  bound  to  H+  mainly  by  polarization 

forces.  A  similar  situation  also  occurs  for  clusters  formed  by  aromatic 

cations  and  rare  gas  atoms.  Since  the  ionization  potentials  of  aromatic 
molecules  are  much  smaller  than  those  of  rare  gas  atoms,  the  positive  charge 
is  located  on  the  aromatic  molecules,  which  thus  attracts  the  rare  gas  atoms 
by  polarization  forces. 

The  brief  survey  presented  above  shows  that  the  study  of  ionic  clusters 
is  a  rap  idly -developing  field  of  research.  Some  of  the  experimental  results, 
especially  the  discovery  of  the  new  (cluster)  class  of  reactions, ^  suggest 
possible  practical  applications  of  ionic  clusters,  such  as  in  catalysis. 

In  this  study  we  shall  present  a  quantum  chemical  treatment  of  the 
relatively  simple  ionic  clusters  R*,  (R^H)*  and  (R^^)*.  These  systems  are 
bound  together  not  only  by  electrostatic  forces  but  sometimes  also  by  valence 
forces  resulting  from  the  electron  (charge)  delocalization,  as  discussed  above 

*4* 

for  R^  ions.  In  order  to  study  theoretically  the  systems  with  ions  where  both 
electrostatic  polarization  and  charge  delocalization  are  of  importance,  we 
have  previously  developed  the  semiempirical  diatomics- in- ionic -systems  (DIIS) 
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9 

method.  This  method  was  first  applied  to  rare  gas -halogen  systems  with 

excited  ionic  states  like  Xe+Ci  and  Xe+HCi",^  and  we  shall  now  use  it  to 

n  n 

calculate  the  electronic  structure  of  rare  gas  and  rare  gas -hydrogen  ionic 
clusters . 

In  the  next  section  we  shall  describe  briefly  the  DIIS  method,  as  well 

as  some  modifications  which  one  needs  to  apply  it  to  the  systems  under 

consideration.  Sections  III  and  IV  deal  with  calculations  of  Ar+,  (Ar  H)+  and 

n  n 

(Ar  H_)  +  and  of  Xe+,  (Xe  H)+  and  (Xe  H„)+,  respectively.  Section  V  given  the 
n  4  nn  n  z 

Conclusions,  and  the  Appendix  presents  the  diatomic  potentials. 


II.  Diatomics - in- Ionic - S vs  terns  (DIIS)  Method 

A.  Closed-Shell  Atoms  with  Electron  Deficiency 

Let  us  consider  a  polyatomic  system  (A^Aj-.-Aj)*  consisting  of  J  closed- 

shell  atoms  or  ions  A^  with  one  electron  deficiency  (hole)  delocalized, 

generally  speaking,  between  these  atoms  (ions).  For  example,  the  neutral 

rare  gas-halogen  system  R^X  is  considered  in  DIIS  as  (R^X ) + .  The  spin  of  a 

system  with  one  electron  deficiency  is  S  -  The  DIIS  wave  function  is 

presented  as  a  linear  combination  of  diabatic  polyatomic  wave  functions  with 

9 

localized  electron  deficiency  as 


M, 


* 


-5  I 

i— 1  m-1 


C. 

im  im 


(1) 


where  is  the  diabatic  wave  function  of  the  configuration  A^A£ .  .  . A^.  .  .Aj , 
with  the  index  m  indicating  the  orientation  of  the  A+  open  shell  (M^  -  3  for  a 
P-symmetry  A*  shell).  For  example,  R^X  is  described  by  nine  diabatic 
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configurations,  namely  R+RX’ ,  RR+X  and  RRX  for  m  -  1,2,3.  The  rare  gas  ionic 
cluster  R*  is  described  by  3n  diabatic  configurations  R...R+...R.  An  ionic 
system  consisting  of  n  rare  gas  atoms  and  one  hydrogen  atom  (R^H)*  can  be  also 
treated  by  the  DIIS  method  if  the  electron  affinity  of  H  is  neglected  and 
singlet  states  only  are  considered.^  The  (R^H)*  system  is  described  by  3n  +  1 
diabatic  configurations,  namely  R. . .R+. . .RH  (m  -  1,2,3)  and  R.  .  .R.  .  .RH+  (m  - 
1). 

The  diagonal  matrix  elements  of  the  DIIS  wave  function  (1)  have  the 
physical  interpretation  as  the  energies  of  the  diabatic  configurations  with 
fixed  localization  of  the  electron  deficiency.  If  we  consider  only  the 
pairwise  interactions  between  atoms  (ions),  then  the  diagonal  matrix  elements 
are  presented  as  a  sum  of  diatomic  interactions, 

5  ’j.j.  *  ui  +  1  Ui-.J  •  J'JrJa-1  •  <2> 

J  j  2  (>J  x)  j 


where  W  is  the  energy  of  the  closed-shell  atom  A.  ,  W.  is  the  diatomic 

J  J  J  2 

potential  between  two  closed-shell  atoms,  A.  -A.  ,  U.  is  the  energy  of  atom  At 

Jl  J2  1  1 

with  an  electron  deficiency,  and  j  is  the  diabatic  potential  between  the 
closed-shell  atom  A^  and  the  atom  with  electron  deficiency,  A*.  The 
transformation  of  the  potentials  between  m-oriented  orbitals  to  the  £  and  II 


potentials  is  described  in  Ref.  9,  where  the  expressions  for  the  off-diagonal 
matrix  elements  are  given  as  well.  The  atom  H  is  considered  here  as  a  closed 
shell  atom.  In  the  case  of  ionic  systems,  like  R+  or  (RnH)+,  the  closed-shell 
atoms  are  neutral,  and  the  atoms  with  electron  deficiency  are  the  ions  R+  or 
H+.  The  diatomic  potentials  between  neutral  atoms,  WDD  and  W_  ,  are  taken  to 
be  known  from  empirical  studies  or  initio  calculations.  The  diabatic 
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potentials  between  neutral  atoms  and  ions,  such  as  R+-R,  R+-H,  R-H+,  have  to 

be  calculated  since  only  adiabatic  potentials,  such  as  W  +  and  W.  +,  are  the 

K-2  CRM) 

physically-defined  values. 

The  diabatic  potentials  ^  can  be  found  by  considering  the  diatomic 
fragments  within  the  DIIS  approximation  and  solving  the  inverse  problem  of  the 
2x2  matrix  eigenvalues,  i.e.,  calculating  the  matrix  elements  for  given 
eigenvalues.  In  the  case  of  the  homonuclear  fragment  R^,  the  2x2  matrix  of 
the  fragment  is 


U  -E  V 

V  U  -  E 


0 


(3) 


where  U  is  the  R+-R  diabatic  potential  and  V  is  the  exchange  term  in  the 

(R+R)-(RR+)  coupling.  Substituting  into  (3)  the  known  adaibatic  potentials 

E(1)  -  22  and  E(2)  -  22  or  E(1)  -  2n  and  E(Z)  -  2n  ,  we  can  easily  find  the 
u  g  g  u’  J 

diabatic  R+R  potential  U  and  the  exchange  term  V  for  2-  or  n-symmetry.  For 
the  case  of  the  heteroatomic  fragment  (RH)  +  ,  there  are  two  different  diabatic 
potentials,  RR+H  and  U^, 


Vh  ‘  E  v 


-  0 


v  URH+  -  E 


(4) 


In  order  to  find  1)^+^,  U^+  and  V  for  two  given  adiabatic  potentials,  for 
example,  -  X^2  and  E^2^  -  B^Z,  one  needs  an  extra  empirical  or  afe  initio 
value  such  as  the  static  or  transition  dipole  moment.  After  the  diabatic 
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potentials  U,  .  and  the  exchange  terns  V.  are  determined  from  Eqs . 

v  im.j  6  ilml’12m2 

(3)  and  (4),  the  diagonal  matrix  elements  (2)  and  the  off-diagonal  matrix 
elements  (presented  in  Ref.  9)  can  be  calculated,  and  the  wave  function  (1) 
and  its  energy  spectrum  can  be  found. 

In  the  wave  function  presentation  (1),  the  spin-orbit  coupling  is  not 
taken  into  account,  at  least  directly.  Indirectly,  the  energy  level  shifts 
resulting  from  the  spin-orbit  coupling  are  incorporated  into  the  system  energy 
via  the  empirical  or  afe  initio  diatomic  potentials  which  are  used  in  the  DIIS 
calculations  as  the  input.  By  neglecting  the  spin-orbit  coupling  in  the  wave 
function,  we  are  losing  some  of  the  excited  states  and  are  excluding  from  our 
consideration  effects  such  as  S-P  mixing  in  collinear  configurations. 

B.  Polarization  Energy 

In  applying  Eqs.  (2) -(4),  we  need  not  care  about  the  polarization 
component  of  the  ion-atom  interactions,  since  it  is  included  in  the  empirical 
or  afe.  initio  potentials  which  we  are  suggesting  to  use.  However,  the 
polarization  of  atoms  gives  rise  to  the  electrostatic  dipole-dipole 
interaction  which  is  of  three-particle  origin  (ion-atom-atom)  and  cannot  be 
considered  in  a  pairwise  way.  The  energy  of  the  dipole -dipole  interaction  is 
usually  small,  although  it  is  not  necessarily  negligible,  and  hence  we  intend 
to  take  it  into  account.  Including  this  energy  as  a  separate  term  A^  within 
the  diagonal  matrix  element  (2)  and  taking  as  zero  the  energy  of  the  separated 
neutral  atoms,  we  obtain 

Him,im“Ii  +  }  \  Wj  J,  +  Y  Uim,j  +  Ai  ’  J  ,J1’ J2  ^  1  ,(5) 

h  J 

where  1^  is  the  ionization  potential  of  the  i-th  atom,  and 
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t-2  I 


jx  jjOj^ 


A.  .  4 

1  •  j  ^  •  j  2 


(6) 


with  A.  .as  the  dipole-dipole  interaction  between  two  neutral  atoms 
*  »  J 1  *  J  o 

14 

polarized  by  the  ion  with  index  i.  Since  the  terms  A.  .  depend  on  the 

1  •  J  ^  >  J  2 

diatomic  polarization  energies  ,  we  must  calculate  them.  P„  are  used  for 
the  A^  calculation  only,  so  that  their  effect  is  supposedly  weak,  which 
reduces  significantly  the  requirements  for  their  accuracy. 

In  the  dipole  approximation  the  ion-atom  polarization  energy  decreases 
with  the  fourth  power  of  the  separation  R, 


Py  “  -Cj/Ry  ,  cj  -  e2a./2  .  (7) 

where  is  the  polarizability  of  the  jrth  atom.  The  dipole  approximation  (7) 
is  justified  provided  that  atoms  do  not  overlap  one  with  another.  From 
physical  reasons  it  is  clear  that  the  overlap  between  an  ion  and  atom  has  a 
damping  effect  on  the  polarization  enefgy  P^ .  We  suggest  taking  this  account 
into  effect  by  multiplying  the  polarizability  in  Eq .  (7)  by  a  damping 
function  y(Ry), 

(Rjj  )  “  aj7(Rj,j)  (8) 


A  suggested  form  for  the  functional  dependence  on  R^  is 


rr«+rn  12-1*1/3 


(9) 


w 


where  and  r^  are  the  ionic  and  atomic  radii,  respectively.  For  the  case  of 

the  interaction  between  polarized  neutral  atoms,  both  and  r^  are  atomic 

radii.  When  the  interatomic  separation  is  larger  than  the  sum  of  atomic 

radii  (non-overlapping  atoms),  the  function  7  is  close  to  unity,  i.e.,  the 

dipole  approximation  is  valid.  When  the  interatomic  separation  is  smaller 

than  the  sum  of  atomic  radii  (overlapping  atoms),  the  damping  function 

decreases  as  ,  so  that  the  ion-atom  polarization  energy  (7)  is  saturated 

4 

with  a  limiting  value  of  -Cj/(r^+r^)  .  The  suggested  form  (9)  for  the  R^ 
dependence  is  of  course  arbitrary,  but  at  least  it  describes  in  a  proper  way 
the  asymptotic  behavior.  If  the  asymptotic  conditions  are  fulfilled,  the  form 
of  the  dependence  of  7  on  does  not  imply  much  about  the  results. 

According  to  our  estimations,  different  reasonable  models  of  the  damping 
function  y(R^)  vary  the  energy  within  a  range  of  0.03  eV  only,  mostly  even 
much  less.  For  the  values  r^  and  r^  in  (9),  we  shall  use  the  atomic  VDW 
radii. 


C.  (R^H^)  System 

The  ionic  rare  gas-hydrogen  system  (R^l^)*  is  described  in  the  DIIS 
method  as  n  closed- shell  atoms  (R)  and  two  open- shell  atoms  (H)  with  one 
positively-charged  hole.  Each  of  the  two  diabatic  configurations  R.  .  .RH+H 
with  a  hole  located  on  one  of  the  H  atoms  is  described  by  a  wave  function  with 
a  fixed  s  -  ^  spin  on  the  only  open-shell  atom  (H)  .  The  situation  is  much 
more  complicated  in  the  diabatic  configuration  R...R+...RHH  when  the  hole  is 
located  on  one  of  the  n  rare  gas  atoms.  In  this  configuration  there  are  three 
open-shell  atoms  (R+,H,H)  with  spins  s  -  which  make  two  different  spin 
states  for  a  total  spin  S  -  Consequently,  each  R.  .  ,R+.  .  .RHH  configuration 
with  fixed  hole  location  is  described  by  six  diabatic  wave  functions  which 
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differ  from  anocher  not  only  by  p-orbital  orientation,  like  in  expansion  (1), 
but  also  by  spin  eigenfunctions.  Taking  into  account  the  hole  delocalization 
between  all  rare  gas  and  hydrogen  atoms,  one  expresses  the  DIIS  wave  function 
as  the  superposition  of  6n+2  diabatic  wave  functions, 

n  3  2  2 

*-}  }  }  Ci»k  *1*  +  }  Cn+J  Vj  •  <10> 

i-1  m-1  k-1  j-1 

where  ^n+j  describes  the  diabatic  configuration  with  the  charge  on  an  H  atom, 

R. . ,RH+H,  9.  .  describes  the  diabatic  configuration  R......R......R,  .  HH  (the 

lmk  °  (1)  (i)  (n) 

numbers  in  parentheses  are  the  atomic  indices)  with  the  m-orientation  of  the 
R*^  P-shell,  and  k  indicates  the  spin  configuration  of  the  three  (R^j  ,H,H) 
uncoupled  electrons.  Following  the  usual  DIM  presentation  of  a  triatomic  wave 
function, ^ we  express  the  diabatic  wave  functions  in  terms  of  one-electron 
orbitals  as 


n 


^ill-yt(lpixPiyPiyPizPizV2l  *  lpixpiypiypizpizsls2 1 } (U) 


*il2  "  76  ^pixpiypiypizpizsls2  ^  ^ pixpiypiypizpizsls2 ^ 


n 


pixpiypiypizpizsls2l)j(;i)^j 


(12) 


j  "  'pjxpjxpjypjypjzpjz'  ’ 


,  j  -  1,2. . .n 


(13) 


n 


Vl  -S2^j  -  *n+2  “  S1  j  ' 


(14) 
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where  Che  bars  denoce  Che  negacive  spin  orlencacion,  p.^,  ply’  piz  are  the  P* 

orbiCals  of  Che  i-ch  rare  gas  aCom,  and  and  Sj  are  Che  Is  orbitals  of  the 

hydrogen  atoms.  The  index  a  -  1  in  *in  anc*  *H2  *n<**-caCes  c^e  x-orientation 

of  the  R.t . .  shell.  The  diabetic  functions  $.  .  for  m  -  2  (y-orientation)  and 
(i)  lmk 

m  -  3  (z-orientation)  are  expressed  by  equations  similar  to  Eqs.  (11) -(12). 

The  p-orbitals  of  the  diabatic  functions  (11) -(13)  are  oriented  along 

fixed  x,y,z-axes.  When  the  diabatic  functions  are  used  to  form  the  H-matrix, 

the  x,y,z  oriented  p-orbitals  have  to  be  transformed  to  2  and  II  orientations 

q 

of  diatomic  fragments.  This  transformation  for  arbitrary  located  atoms  is 
performed  by  an  S  matrix  whose  elements  are 


ij.xS  “ 

<Y 

•x.)/R..  , 

sij  ,xn  ~ 

riJ/Rij  '  Sij  ,x»  "  0  1 

,  (15a) 

ij.ys  " 

<yj- 

-V/Rij  ■ 

sij.yn  " 

• 

sij»yx  ” 

<rj'zl)/riJ  ’ 

(15b) 

ij.zZ  ’ 

<2j 

'Zi)/Rij  ’ 

sij,zn  " 

(Zj  * z^) /(R^j  ) 

l 

Sij  ,xx  " 

(15c) 

2  2 

rij  ■  +  <zj’ ^ 


(16) 


where  x^ ,  y^  and  z^  are  the  i-th  atom  coordinates,  R^j  is  the  interatomic 
separation,  and  the  H-orientation  of  the  R*^  P- shell  is  indicated  by  II  when 
located  in  the  plane  formed  by  and  the  x-axis,  and  by  ir  when  perpendicular 
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to  this  plane.  We  shall  also  use  later  the  following  combinations  of  the  S 
values : 


T 


ij  ,m 


Sij,mn 


+  S 


2 

ij  ,mjr 


(17a) 


“  Sijimin  Sij,m2n  +  Sij,mlJr  Sij,m2w  *  (i7b) 

9  47  58 

Applying  both  the  DIIS  and  DIM  ’  approaches  and  using  the  diabatic 

functions  (11)-(14),  we  obtain  the  (6n+2)x(6n+2)  H*matrix  which  we  shall 

simply  display  without  going  into  the  details  of  its  construction.  The 

diagonal  matrix  element  for  R .. . . . .R^, .  . . .R .  .H.  is 

(1)  (i)  (n)  (n+1)  (n+2) 


n-1  n 


n 


Himk,imk  “  lK  +  )  )  Wj1j2+I  [Sij,mE  Uij,2  +  Tij,m  Uij,n] 

1  /  J  V  J  \  L.  L.  •  1 


jl  <W  ‘  “  J“l 

n  2  2 

^  \  \  +  }  fSi,n+j  , 


+ 

j'x-i  j-i  ‘  j-i 

+  Wn+l,n+2  ’  J » J i » J 2  *  1  *  k  “  1'2  - 


i , n+ j , Z  4  i,n+j,Z  4  i,n+j,Z 


0(k)  .  T  0(k)  , 

m2  i,n+j,Z  i,n+j,m  i,n+j,nJ 


0<2>  -  a  h 

i,n+j,Z  4  u i , n+  j ,  Z  r  4  wi,n+j,Z  ’ 


+  f  V 


(18) 


(19a) 


(19b) 


(k) 

and  the  same  for  n+j  jj.  In  Eqs.  (18)-(19),  1^  is  the  rare  gas  ionization 

potential,  W.  .  are  the  R-R  potentials,  W.  are  the  R-H  potentials, 

J  ]_J  2  ^  1  ’ 
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Wn+1  n+1  are  C^e  H’H  (k  “  1)  and  32  (k  -  2)  potentials,  ^  and  ^ 

are  the  R*  -R..  diabetic  potentials,  and  ^U.  _  „  (  u.  .  _)  and  3U.  .  _ 

(i)  (j)  i,n+j,2  i,n+j,ir  i,n+j,Z 

(  UL  n+j  are  the  Z- symmetry  (II-symmetry)  R+-H  diabatic  potentials  in  the 
singlet  and  triplet  states,  respectively.  The  diagonal  matrix  element  for 

R> • •RH(n+l)H(n+2)  is 


n-1  n 


H 


n+1 , n+1 


-:H  +  )  I  \h*l  <Hj. 
J1  <J2(>J1) 


n+2  +Un+l,j>  +Un+l,n+2  ’(20) 


and  the  same  for  H  ,  In  Eq .  C20),  U  ,  .  are  the  H  -R  diabatic 

n+3 , n+2  n  n+1 , j 


potentials,  and  Un+^  n+2  Is  the  R(n+2)  d^a^at^-c  potential 

The  off-diagonal  matrix  element  between  two  R^  ,  . .R+^  •  • 
configurations  with  different  (k  -  1,2)  spin  orientations  are: 


H  _  H 

iml,im2  4 


j-1 


3« _ -  V  ...  „>I 


“i.ntj  ,-r) 


+  T  ( 

i,n+j  ,m  “i,n+j,n  ''i.n+j.II 


(21) 


The  off-diagonal  matrix  elements  between  two 

R^j . . -R^jj  •  •  'R(n)H(n+2jH^n+2^  configurations  with  different  (m^,m2)  orbital 
orientations  are 


n 


irn^k, im2k 


'  )  ^iJ.m.S  Sijm2S  Uij,E  +  Tij  ,m  m,Uij  ,II] 
j-1  L  Z  1  1 
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c 

\  r s  g  0Ck>  ,  -r>  0(k)  i  1  »•  1 

i.n+j.n^Z  ai,n+j,m22  uij  ,S  +  i,n+j  ,nJ  ’  J  L 


(22) 


c 

Him^l , im22  , im^l  4^  ^  ^  ^i,n+j ,  m^Z  S i , n+ j , m^ 


X  ( -  ^TJ  \  4.  t  ( ^TT 

i . n+ j  , 2  i,n+j,Z;  i.n+j.n^m^  i,n+j,n 

•  Vi,..11  • 


(23) 


The  off-diagonal  matrix  elements  between  the  configurations 

R(l) • • •R(i1) • * 'RnH(n+l)H(n+2)  and  R(l) *  * *R(i  ) 1 1 -RnH(n+l)H(n+2)  are 


i^®^l  >  ■i'2n,2^’  Hnl^’  *2°2^  i^i2S 


+  ^Sili2ml11  Si2ilm211  +  Sili2ml,r  Si2iim21r^Vtli2,r 


Hil0ll,i2o22  “  0  ’  (25) 

where  V  and  V  .  are  the  exchange  terms  of  the 

1lt2i'  1  2a 


(R(i  )^*^R(i  )R|i  )^  coupling  (see  Eq.  (3)).  The  off-diagonal  matrix 

elements  between  the  configuration  R,  .H,  -  *  and  the 

(1)  (i)  (n)  (n+1)  (n+1 ) 

configurations  and 

R(l) ’ * ,R(i) ■ ‘ -R(n)H(n+l)H(n+2)  are 


Himl , n+j  “  " ^ i , n+ j , m , Z  Vi,n+j,2  ’  j  “  1,2 


(26) 
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Himl,n+j  •  J  *  l'2 


(27) 


where  V.  .  _  is  the  exchange  term  in  the  (R+H)-(RH+)  coupling  (see  Eq.  (4)) 
i ,  n+J  ,  2. 


The  n-component  of  the  coupling  is  zero 

.+ 


The  off-diagonal  matrix  element 

+ 


between  the  configurations  *(1) • •  R(n)H(n+L)H(n+1)  and  »(1)-R(n)H(nfl)H(n|2) 
is  equal  to  the  exchange  term  of  the  (H+H)-(HH+)  coupling  in  the  molecule 
(see  Eq.  (3)), 


^n+1 , n+2  ^n+1 , n+2 


(28) 


When  the  delocalization  of  the  charge  is  neglected  (i.e.,  setting  the  exchange 
terms  V  equal  to  zero),  we  obtain  for  R...R+...RHH  the  usual  DIM  matrix 
elements . 

The  diabatic  potentials  0. ,  _  and  U. .  of  the  R+R  diatomic  fragments 

1J  ,2i  ij 

(Eqs.  (18)  and  (22)),  as  well  as  n+2  of  the  H+H  fragment  (Eq.  (20)),  are 
found  by  the  inverse  solution  of  Eq.  (3),  as  described  in  Subsection  II. A 
above.  The  situation  is  much  more  complicated  for  the  diabatic  potentials 

\»),v  <23»  *”d  Vi.j  <E<<-  (20))  of  the 

R  H  and  RH  diatomic  fragments.  In  our  approach,  the  (RH)  diatomic  fragment 
is  described  by  a  3  x  3  matrix  which  is  obtained  by  removing  to  infinity  all 
other  atoms.  Using  Eqs.  (18),  (20),  (26)  and  (27),  we  have 


1  Itt  .2  3u  .  E 

4  “i,n+l  4  Ui,n+1  E 


*^(^U  -*U  ) 

4'  i,n+l  i,n+l' 


•^(^U  -^TJ 
4'  i,n+l  ui, 


n+l'  4  u  i ,  n+1  4  ui,n+l 


i,n+l 


-yi  V 

1 1 


n+1 


i,n+l 


-73  V 


i  ,n+l 


U 


n+1 , i  -  E 


-  0  , 


(29) 
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where  we  have  simplified  the  notation  by  dropping  the  subscript  Z  for  the 

fragment.  Using  the  (RH)+  potentials  and 

(asymptotic  to  R*+H)  and  ^  (asymptotic  to  R+H*)  as  the  eigenvalues 

(2)  (3^ 

E  and  E  for  Eq .  (29),  we  can  find  numerically  the  diabetic  potentials 

*U.  ...  ^U.  -  and  U  ,  The  construction  of  the  exchange  term  V,  .  of 

i,n+l’  i,n+l  n+l,i  i,n+l 

the  (RH)*  diatomic  fragment  was  described  in  Subsection  II. A  (see  Eq.  (4)). 

The  diabatic  n-potentials  Hj.  .  _  and  ^U.  .  _  coincide  with  the  adiabatic 

*  i,n+j,n  i.n+j.H 

1  3 

W,  .  _  and  W,  .  .  _  potentials,  since  the  exchange  term  V.  ,  _  is  equal 
i,n+j,n  i ,  n+ j ,  II  r  6  r.n+j.H  H 

to  zero.  The  adiabatic  diatomic  potentials  of  the  diatomic  fragments  R^ ,  > 

RH,  R*,  and  (RH)*  for  R  ■  Ar,Xe  are  presented  in  the  Appendix. 


III.  Argon  and  Argon- Hydrogen  Ionic  Clusters  Ar*,  (Ar^H)*  and  (Ar^H^)* 

Considering  the  structure  of  rare  gas-hydrogen  ionic  clusters,  we  have 

decided  to  look  at  Ar  and  Xe,  since  they  are  expected  to  form  different  kinds 

of  clusters  with  H  due  to  the  difference  in  their  ionization  potentials. 

Specifically,  the  Ar  ionization  potential  (IAr  ”  15.76  eV)  is  much  higher  than 

the  H  ionization  potential  (I„  -  13.6  eV) ,  whereas  the  Xe  ionization  potential 

(I^e  -  12.13  eV)  is  smaller  than  that  of  H,  which  leads  to  different  charge 

distributions  in  the  argon-hydrogen  and  xenon -hydro gen  clusters.  But  before 

considering  the  rare  gas-hydrogen  systems,  we  shall  present  the  results  of 

calculations  on  the  rare  gas  ionic  clusters  Ar+  and  Xe*. 

°  n  n 


A.  Ar*  Clusters 
_ Q _ 

The  results  of  the  calculations  of  the  small  ionic  clusters  Ar^ , 
Ar£,  Ar*  and  Ar^  are  presented  in  Table  1  and,  partly,  in  Fig.  1.  We  shall 
now  compare  these  results  with  the  known  experimental  and  theoretical  results 
which  have  been  obtained  for  the  Ar*  clusters. 
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4* 

According  Co  experimencal  studies  of  Che  simplest  cluster,  Ar^,  its 

energy  of  dissociation,  Ar+  -*  Ar+  +  Ar,  is  about  0.2  eV.  ^ ^  ^  ^  There 

are  two  alternative  structures  of  the  Ar^  cluster,  namely  an  asymmetrical 

Ar+Ar  structure  with  the  charge  concentrated  on  Ar^  and  a  structure  with  the 

charge  distributed  among  all  three  atoms.  In  the  case  of  the  Ar^Ax  structure, 

where  the  neutral  atom  is  bound  to  the  Ar^  molecule  by  polarization  forces, 

the  triangle  geometry  is  expected  to  be  the  most  stable  one.^  When  the 

charge  is  delocalized  among  all  three  atoms,  the  exchange  interaction 

contributes  significantly  to  the  binding  so  that  the  symmetrical  linear 

geometry  is  expected  to  be  the  most  stable  configuration.  The  first  quantum 

chemistry  calculation,  performed  by  the  approximate  method,  gave  preference 

6  2 

to  the  asymmetrical  triangle  structure.  However,  the  Cl  S^L  initio 

24  25  27  29 

calculations,  ’  as  well  as  semiempirical  calculations,  ’  found  the 

symmetrical  linear  configurtaion  to  be  the  most  stable  one.  The  experimental 

spectral  data,  unfortunately,  cannot  provide  direct  evidence  about  the  Ar^ 

geometry.  The  conclusions  made  by  the  indirect  analysis  of  the  experimental 

data  are  contradictory.  Whereas  some  of  the  experimental  works  confirm  the 

symmetrical  linear  Ar^  structure , ^ ^ other  napers  substantiate  the 

asymmetrical  Ar^Ar  structure,  either  of  triangular^  or  linear^ ’ ^  geometry. 

Such  discrepancy  in  the  conclusions  is  due,  most  probably,  from  the 

features  of  the  Ar*  potential  energy  surface.  Since  Ar^  has  a  very  flat  Ar£- 
19  28  64 

Ar  potential,  ’  ’  the  vibrational  motion  shifts,  for  the  most  part,  the 

Ar^  configuration  far  away  from  the  equilibrium  configuration,  thus 

significantly  affecting  the  electronic  transitions.  This  problem  has  been 

28 

considered  in  detail  by  Gidea  and  Amarouche,  who  calculated  the  potential 
energy  surfaces  of  Ar^  by  a  method  similar  to  DIIS.  A  trajectory  study 
performed  on  these  potential  energy  surfaces  shows  that  the  symmetrical  linear 
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Ar^  minimum  energy  configuration  alone  provides  the  spectrum  close  to  the 

28 

experimental  one,  if  one  takes  into  account  the  vibrational  motion. 

In  our  calculations  the  triatomic  cluster  Ar^  has  a  symmetrical  linear 

24  -  28 

equilibrium  configuration,  in  accord  with  other  calculations.  The 

2.-*  2  transition  energy  in  the  equilibrium  configuration  is  found  in  our 

Kg 

calculations  to  be  equal  to  2.26  eV,  close  to  the  other  theoretical  values  of 
24  28 

2.36  eV  and  2.34  eV.  Our  2.  -*  2  transition  moment  is  8.8  D,  compared  to 

K  S 
24 

the  fib.  initio  value  of  8.24  D.  For  the  Ar-Ar  distance  at  the  equilibrium 

configuration,  our  calculation  gives  R  -  2.59  A,  exactly  the  same  as  the 

28 

semiempirical  calculation  of  Gadea  and  Amarouche  but  slightly  less  than 

25  4*  4- 

initio  calculations  (2.62  A).  The  Ar^  -+  A^  +  Ar  dissociation  energy  in  our 

calculation  is  D  —  0.203  eV,  close  to  the  experimental  values  of  0.22  eV^’^ 

21 

and  0.18  ±  0.05  eV.  Other  calculations  give  values  a  little  bit  lower  than 

24  25  27 

our  dissociation  energies,  namely,  0.18  eV,  0.16  eV,  0.20  eV  and  0.17 


The  most  stable  Ar*  cluster,  according  to  our  calculations,  is  formed  by 

an  Ar*  ion  and  almost  neutral  (q  -  +0.002)  Ar  atom  separated  from  the  nearest 

charged  atom  by  3.68  A  (see  Table  1  and  Fig.  l.III).  The  dissociation  energy 

for  the  Ar*  -*  Ar*  +  Ar  detachment  process  is  found  to  be  equal  to  D  -  0.047 

eV,  a  typical  energy  of  the  polarization  attraction.  Another  Ar*  cluster  is 

found  in  an  asymmetrical  linear  configuration  with  one  almost  neutral  (q  - 

+0.01)  atom  separated  from  the  nearest  atom  of  Ar^  by  3.3  A  (see  Table  I  and 

Fig.  l.V).  The  dissociation  energy  of  this  quasistable  configuration  is  D  - 

0.031  eV,  and  the  potential  barrier  of  the  transition  to  the  most  stable 

27 

configuration  (Fig.  l.III)  is  about  0.001  eV  only.  In  the  DIM  calculation,  v 
the  linear  configuration  with  the  dissociation  energy  of  0.043  eV  is  found  to 
be  the  most  stable  one.  According  to  experimental  studies, the  spectrum 
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of  Ar+  and  other  small  Ar+  clusters  is  similar  to  the  Ar^  spectrum,  so  that 

the  conclusion  is  made  that  the  Ar+  (n  >  3)  clusters  have  a  structure  of 

n 

Ar^Arn  This  experimental  finding  supports  our  result  that  Ar^Ar  is  the 

most  stable  structure  for  the  Ar+  cluster.  A  similar  result  is  obtained  for 

4 

Ar ^  and  Ar+  (see  below).  Like  in  the  experimental  studies,  the  spectrum  of 

the  most  stable  Ar+  cluster  is  found  to  be  similar  to  that  of  Ar^  (Table 

2. III).  It  consists  of  two  transitions  with  large  transition  moments  and 

2  + 

transition  energies  close  to  one  another  and  to  the  2  excitation  energy  of 

S 

Ar^.  so  that  they  give  one  absorption  band  with  its  center  close  to  that  of 
the  Ar^  absorption.  It  is  interesting  to  note  that  in  the  Ar*  excitation 
states  a  strong  charge  transfer  takes  place,  which  makes  it  impossible  to 

Save 


describe  the  excited  (Ar^)  cluster  as  an  (Ar^)  Ar  structure  (Table  2) 

2  +  2  +  + 

for  the  transitions  which  resemble  the  2  ■+  2  transition  of  Ar, ,  we  find 

u  g  3’ 

also  in  Ar*  a  low-energy  (1.79  eV)  transition  with  a  small  transition  moment 
of  /t  -  0.7  D  (Table  2). 


In  addition  to  the  Ar^  clusters  with  the  ground- state  Ar^Ar  structure, 
we  find  a  quasistable  symmetrical  configuration  with  all  atoms  bearing  the 
same  charge  of  q  -  +0.25,  (Table  l.IV).  In  this  configuration,  Ar£  has  the 
geometry  of  a  regular  pyramid  with  a  separation  of  2.836  A  between  every  two 
atoms.  The  energy  of  this  fully- symmetrical  structure  is  only  0.1  eV  above 
the  energy  of  the  most  stable  Ar*  cluster.  The  spectrum  of  the  symmetrical 
Ar£  cluster  is  quite  different  from  that  of  Ar^.  The  transition  energy  to  the 
three  lowest  excited  states  which  form  a  degenerate  level  is  only  1.46  eV, 
compared  to  2.26  eV  for  the  first  allowed  transition  in  Ar^  (Table  2).  The 
symmetrical  Ar£  configuration  is  separated  from  the  stable  Ar^Ar  configuration 
by  a  high  barrier  of  roughly  0.25  eV. 
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4*  4* 

The  Ar,.  and  Arg  clusters  are  formed  in  a  similar  way  as  the  stable  Ar^ 
cluster,  i.e.,  by  a  Ar^  ion  and  neutral  or  almost  neutral  Ar  atoms.  The 
energy  of  one  Ar  atom  detachment  is  0.05  eV  for  Ar+  and  0.06  for  Ar+  (Table 
1).  The  Ar^  and  Ar£  spectra  are  similar  to  the  spectrum  of  the  stable  Ar* 
cluster  (Table  2.  III).  The  experimental  studies  confirm  that  Ar*  is  an  ionic 
core  in  small  Ar*  clusters .  ^ ^  However,  in  large  Ar*  clusters,  n  >  15,  the 
Ar^  molecule  was  found  to  be  the  most  stable  ionic  core.^ The  transition 
from  the  Ar^  to  Ar^  core  with  an  increase  of  cluster  size  can  be  explained  by 
the  polarization  effects.  The  energy  of  the  neutral  atoms  polarization  by  the 

+  4- 

ionic  is  slightly  larger  for  the  Ar^  than  Ar^  core  because  of  the  difference 
in  their  dimension.  In  small  Ar*  clusters,  this  difference  in  polarization 
energy  is  small  compared  to  the  Ar^  dissociation  energy,  and  consequently  it 
is  not  expected  to  affect  the  core  structure.  In  large  Ar*  clusters,  the 
polarization  energy  is  close  to  that  in  solid  Ar  where,  according  to  our 
estimations,  the  difference  between  the  polarization  energies  of  the  Ar^  and 
Ar^  cores  lies  somewhere  between  0.3  and  0.6  eV.  This  energy  difference  is 
larger  than  the  Ar^  dissociation  energy,  so  that  the  Ar^  core  may  become  more 
stable  than  the  Ar^  one. 

B.  (Ar  H)+  Clusters 
_ Q _ 

The  structure  of  the  (ArH)*  clusters  is  simple:  all  of  them  are  formed 
by  an  (ArH)*  ionic  molecule  and  n-1  neutral  Ar  atoms  which  are  bound  to  (ArH)* 
by  polarization  forces  (Table  3).  Since  in  the  (ArH)+  molecule  the  H  atom 
bears  more  charge  than  the  Ar  atom  (q^  -  +0.569),  the  neutral  Ar  atoms  are 
located  on  the  H  side  of  the  (ArH)+  ion  (Fig.  2) .  The  simplest  (ArnH)* 
cluster,  Ar(ArH)*,  has  a  bent  geometry  with  the  neutral  Ar  atom  at  a  distance 
of  2.81  A  from  H.  In  the  ArjCArH)*,  Ar^(ArH)*  and  Ar^(ArH)*  clusters,  the 
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neutral  atoms  are  located  symmetrically  around  the  (ArH)+  molecular  axis  at 
the  distances  of  2.79  A,  2.79  A  and  2.84  A,  respectively.  In  the  Arg(ArH)+ 
clusters,  two  of  the  Ar  atoms  are  more  tightly  bound  to  (ArH)+  than  the  other 
four  Ar  atoms.  In  the  sequence  of  Ar^(ArH)+  clusters,  the  H  charge  increases 
monotonically  with  k.  The  energy  of  an  Ar  atom  detachment  from  an  Ar^(ArH)  + 
cluster  is  relatively  large,  varying  from  0.09  eV  for  k  -  6  to  0.19  eV  for  k  - 
3.  The  electronic  spectrum  of  the  Ar^(ArH)+  clusters  is  practically  the  same 
as  the  spectrum  of  the  isolated  (ArH)+  molecule.  The  Ar^(ArH)+  clusters  can 
probably  be  detected  by  studying  their  rotational  spectrum,  since  the 
Ar^(ArH)+  dipole  moment  depends  strongly  on  k  (Table  3) . 


C.  (Ar  H_)+  Clusters 
n  Z _ 

The  ionization  potential  of  Ar  (15.76  eV)  is  close  to  the  ionization 

potential  of  the  molecule  (15.426  eV)  ,  which  leads  to  strong  coupling 

between  the  states  Ar+H2  and  ArH^.  This  coupling  affects  significantly  the 

potential  energy  surfaces  of  the  (ArH2)+  system,  which  have  been  calculated  by 
47 

the  DIM  method.  We  shall  not  consider  here  the  (ArH2>  potential  energy 
surfaces,  i.e.,  we  shall  not  consider  the  energy  E  as  a  function  of  the 


interatomic  distances  R 


Ar-H, 


Ar-H . 


and  R 


H-H’ 


but  restrict  our  task  to 


the  study  of  ionic  clusters  formed  by  the  Ar  atom  or  atoms  and  the  H2 
molecule.  The  calculation  of  the  (ArnH2)+  system  is  performed  by  the  modified 


DIIS  method  described  in  Subsection  II. C  above. 


The  (ArH2)+  clusters  can  be  formed  by  Ar+  +  H2  (-4.746  eV)  ,  Ar  +  H2 
(-4.95  eV)  and  H  +  (ArH)+  (-6.215  eV) .  It  follows  that  the  most  stable 
(ArH2>+  cluster  is  expected  to  be  formed  by  the  (ArH)+  molecule  and  H  atom. 
However,  the  polarization  attraction  of  the  H  atom  to  the  (ArH)+  ion  is  so 
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weak  that,  according  to  our  calculation,  no  (Arl^)*  cluster  can  be  formed. 

The  ionic  (ArHj)  cluster  can  be  formed  in  a  quasistable  configurations  only. 

We  find  that  the  quasistable  (ArHj)*  cluster  with  the  lowest  energy  is  formed 

by  an  almost  neutral  Ar  atom  (q  -  0.02)  and  the  molecular  ion  with  an  H-H 

separation  of  1.05  A,  slightly  shorter  than  in  the  isolated  (1.06  A).  The 

Ar  atom  is  located  in  a  plane  perpendicular  to  the  H*  axis  at  a  distance  of 

2.935  A  from  the  H+  center  (Fig.  2).  The  energy  of  the  Ar  detachment  is 

relatively  large,  D  -  0.193  eV.  The  ArH^  cluster  is  separated  from  the  stable 

state  of  separated  (ArH)+  and  H  by  a  high  barrier  of,  roughly,  0.75  eV.  The 

saddle  point  for  the  ArH*  ■*  H  +  (ArH)+  dissociation  lies  in  the  collinear 

configuration  with  the  interatomic  distances  ^  =  1.5  A  and  »  1.2  A. 

47 

Unfortunately,  in  the  DIM  calculation  the  geometry  of  Fig.  2  was  not  studied 

so  we  do  not  have  a  comparison  for  our  results. 

The  Ar  H*  clusters  with  a  few  Ar  atoms  have  a  structure  similar  to  that 
n  z 

of  Arf^  (see  Fig.  2  and  Table  4).  In  all  of  these  clusters  up  to  n  -  5 ,  the 
Ar  atoms  are  located  in  a  plane  perpendicular  to  the  H+  axis  at  distances  of 
about  3  A.  The  energy  of  the  Ar  detachment  is  0.19-0.2  eV,  slowly  decreasing 
with  n  for  n  >  2.  The  common  charge  of  the  Ar  atoms  slightly  increases  with  n 
reaching  q  -  ±0.075  in  Ar^H^.  The  decrease  of  the  charge  with  n  leads  to 
some  decrease  of  the  H-H  interatomic  distance,  which  is  1.03  A  in  Ar^H*.  The 
larger  Ar^^  clusters  (n  >  5)  are  expected  to  be  much  less  stable  than  the 
clusters  presented  in  Table  4.  The  photoabsorption  spectrum  of  the  Ar^H^ 
clusters  is  very  different  from  the  H^  spectrum,  since  the  lowest  excited 
states  are  built  up  by  charge  transfer  to  Ar  atoms  and  the  formation  of  Ar^Hj 


(Table  4). 
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+  +  + 

IV.  Xenon  and  Xenon -Hydro gen  Clusters  Xe^,  (Xe^H)  and  (Xe^H^) 

A.  Xe+  Clusters 

n _ 

The  results  of  the  calculations  of  the  small  ionic  clusters  Xe^,  Xe^, 
Xe^  and  Xe^  are  presented  in  Table  5  and  (partly)  in  Fig.  3.  In  our 
calculations,  the  triatomic  cluster  Xe^  has  a  symmetrical  linear  structure 
with  50%  of  the  charge  at  the  central  atom,  in  accord  with  other 

0  /  O  Q 

calculations.  ’  Our  calculation  gives  3.38  A  for  the  Xe-Xe  separation  in 

Xe^.  compared  to  3.47  A  and  3.32  A  in  other  calculations.  The  Xe^  -*■  + 

Xe  dissociation  energy  in  our  calculation  is  D  -  0.20  eV,  which  is  0.07  eV 

22 

lower  than  the  experimental  value  of  0.27  eV.  Our  dissociation  energy  is 

24  26 

almost  exactly  halfway  between  the  values  of  0.12  eV  and  0.36  eV  obtained 

2  +  2  + 

in  two  other  calculations.  For  the  2  -*■  2  transition  energy  from  the 

u  g 

equilibrium  point,  we  obtain  T  -  1.87  eV,  a  little  bit  higher  value  than  in 
Ref.  24  (1.60  eV) . 

Our  results  obtained  for  Xe+  clusters  with  n  >  3  are  quite  different 

n 

from  those  obtained  in  Ref.  29.  According  to  our  calculations,  the  most 
stable  Xe^  cluster  has  a  structure  of  a  regular  pyramid,  like  Ar^  (Fig.  l.IV), 
and  the  lowest  energy  quasistable  Xe^  cluster  has  a  structure  Xe^Xe  with  the 
weakly -charged  atom  (q  -  0.014)  located  to  the  side  of  the  Xe^  axis  (Table  5, 
Fig.  3).  We  find  also  a  quasistable  linear  symmetrical  Xe^  cluster  which 
energy  lies  0.012  eV  above  the  Xe^Xe .  In  contrast  to  this,  in  Ref.  29  the 
most  stable  Xe^  structure  (Xe  detachment  energy  of  0.159  eV)  Is  just  the 
linear  symmetrical  one.  A  structure  like  our  regular  pyramid  (Fig.  3. Ill)  was 
not  found  in  Ref.  29  at  all.  We  find  the  collinear  Xe^  cluster  V  to  be 
separated  from  the  bent  geometry  Xe^Xe  cluster  IV  by  a  small  barrier  of  0.03 
eV.  The  Xe^Xe  cluster  IV  which  can  be  formed  easily  as  the  result  of  the 
attachment  of  Xe  to  Xe^  (attachment  energy  of  0.08  eV)  is  separated  by  a  high 
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barrier  of  approximately  0.3  eV  from  the  most  stable  Xe*  cluster,  III.  It 
follows  that  the  Xe^X  cluster  can  exist  for  a  relatively  long  time  when  the 
temperature  is  not  high. 

The  clusters  Xe^+  and  Xe^  are  presented  in  Table  5  for  both  the  Xe^.e^  ^ 
and  Xe^Xen  3  structures.  The  energy  of  Xe  detachment  varies  in  the  interval 
0.07-0.11  eV.  The  binding  energy  of  the  attached  atoms  is  derived  from 
polarization  forces,  since  these  atoms  are  almost  neutral. 


In  the  family  of  XejXe^  3  clusters  (IV,  VII  and  IX  in  Table  5),  the 

2  +  + 

strongest  transition  is  to  the  excited  state  similar  to  2^  of  Xe^,  i.e.,  with 
a  very  small  charge  on  the  central  Xe+  atom.  In  some  of  these  clusters,  an 
important  portion  of  the  charge  is  transfered  in  this  excited  state  from  Xe* 
to  the  attached  atoms,  for  example,  q  -  0.38  in  the  Xe*  cluster  IV.  As  in  the 
case  of  the  Ar*Arn  3  clusters,  the  transition  energy  is  not  affected  much  by 
the  charge  transfer,  so  that  in  the  Xe3Xn3  clusters  a  strong  photoabsorption 
is  expected  in  the  same  range  as  in  Xe3  .  In  addition  to  this  2^  -*■  2^-like 
transition,  more  allowed  transitions  with  smaller  transition  energy  are  found 
in  the  Xe^Xe^  3  clusters  (Table  6). 

In  the  symmetrical  Xe^  configuration  III,  there  is  threefold  degenerate 
level  with  an  energy  of  1.20  eV  above  the  ground  state  (and  a  transition 


moment  of  2.7  D) .  In  Xe^Xe  and  Xe^Xe2 ,  two  of  the  excited  states  remain 
degenerate  with  almost  the  same  excitation  energy  as  in  Xe*,  whereas  the  third 
state  excitation  energy  is  lowered  by  about  0.15  eV.  The  linear  Xe^  cluster 


has  only  one  strong  transition  2^  ■*  2^,  whose  energy  is  1.46  eV,  about  0.4  eV 

less  than  in  Xe3Xe .  The  differences  between  the  spectra  of  Xe*  III  (regular 

pyramid),  Xe3Xe  IV  and  Xe+  V  (linear)  may  help  in  the  experimental  search  for 

the  most  stable  Xet  cluster. 

4 
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B.  (Xe  H)+  Clusters 

As  in  the  case  of  the  (Ar^H)*  clusters,  the  (Xe^H)*  clusters  are  formed 
by  the  (XeH)+  molecule  which  attracts  more  neutral  Xe  atoms.  Since  in  the 
(XeH)+  molecule  the  Xe  atom  bears  most  of  the  charge,  the  neutral  Xe  atoms  are 
located  on  the  rare  gas  atom  side  of  (XeH)+  (Fig.  4).  The  Xe(XeH)+  cluster  is 
found  to  have  a  linear  geometry  with  a  Xe  +  (XeH)  dissociation  energy  of  0.046 
eV,  which  is  much  less  than  the  Ar(ArH)+  dissociation  energy  (0.163  eV)  .  The 
energy  of  the  Xe  detachment  in  the  Xe^(XeH)+  clusters  varies  in  the  interval 
0.04-0.08  eV  (Table  7).  The  charge  of  H  in  (XeH)+  slightly  decreases  when  the 
number  of  attached  neutral  atoms  increases.  The  spectrum  of  the  Xe^(XeH)  + 
clusters  is  practically  the  same  as  that  of  isolated  (XeH)+. 

C.  (Xe  H„)+  Clusters 

n  2 _ 

The  most  stable  (Xe^i^)*  clusters  are  expected  to  be  formed  by  charged 
Xe  atoms  and  a  neutral  molecule,  since  the  ionization  potential  is  about 
3.3  eV  higher  than  that  of  Xe.  The  calculation  of  the  simplest,  (Xel^)*, 
system  shows  the  existence  of  a  strongly-bound  cluster  with  the  energy  of  the 
Xe+  +  H£  dissociation  of  D  -  0.456  eV.  This  cluster  has  a  linear  geometry 
with  the  Xe-H  distance  of  1.79  A,  0.18  A  larger  than  in  the  (XeH)+  molecule, 
and  the  H-H  distance  of  0.81  A,  0.07  A  larger  than  in  molecule  (Fig.  4). 
Almost  a  whole  charge  is  concentrated  on  the  Xe  atom  (q  -  0.955),  but  the 
valence  bond  between  Xe+  and  is  strong.  The  excitation  energies  for  the 
allowed  transitions  (in  the  equilibrium  configuration)  are  10.76  eV  and  18.56 
eV,  whereas  in  Hj  there  is  one  transition,  19.2  eV,  and  in  (XeH)+  there  are 
two  transitions  but  with  energies  of  7.93  eV  and  13.3  eV. 

The  structure  of  (Xe^j)*  clusters  with  more  than  one  Xe  atom  is  quite 
different  from  that  of  (XeHj)*  due  to  the  formation  of  Xe+  or  Xe^  molecules. 
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All  (Xe  H0)+  clusters  that  we  have  studied  (n  >  1)  have  the  structure  Xe+H0 
n  z  n  i. 

with  a  practically  neutral  H2  molecule.  In  the  most  stable  Xe+Hj  cluster,  the 
axis  of  both  molecules  are  perpendicular  to  one  another  (Fig.  4).  The 
dissociation  energy  of  this  cluster  is  0.084  eV.  In  the  Xe+H2  clusters  three 
Xe  atoms  form,  in  contrast  to  the  free  Xe+  cluster,  a  triangular  Xe*  system 
with  Xe-Xe  distances  in  the  interval  3.50-3.58  A.  In  one  of  the  Xe*H2  isomers 
with  the  Xe2H2  -*  Xe*  +  H2  dissociation  energy  of  D  -  0.108  eV ,  the  Xe*  plane 
is  perpendicular  to  the  H2  axis  with  its  center  located  on  this  axis  at  3.08  A 
from  the  H2  center.  In  another  isomer  (D  -  0.097  eV)  ,  H2  is  located  in  the 
Xe*  plane.  In  contrast  to  the  Xek(XeH)*  clusters,  in  Xe*H2  clusters  (n  >  2), 
the  valence  binding  affects  to  a  certain  degree  the  cluster  structure. 

V.  Conclusions 

9 

1.  The  semiempirical  diatomics* in-molecules  (D1IS)  method  proves  to  be 
an  effective  tool  to  treat  rare  gas  R*  and  rare  gas-hydrogen  (R^H)*  ionic 
clusters.  The  DIIS  results  obtained  for  the  triatomic  Ar*  and  Xe*  clusters 
are  close  to  those  obtained  in  other  studies. 

2.  According  to  our  calculation,  there  are  two  very  different  isomers 
of  R*  clusters.  One  has  the  structure  of  a  R*  ion  and  almost  neutral  R  atom, 
whereas  the  other  isomer  has  the  symmetrical  geometry  of  a  regular  pyramid 
with  all  four  atoms  bearing  the  same  charge.  We  find  the  Ar^Ar  and  the 
symmetrical  Xe*  clusters  to  be  the  most  stable  isomers,  so  that  Ar*  (n  -  4-6) 
and  Xe*  (n  -  5,6)  to  have  ionic  cores  Ar^  and  Xe£,  respectively.  The 
experimentally  detected  transition  from  the  Ar*  to  Ar*  ionic  core  in  large 
(n  >  15)  Ar*  clusters  results,  most  probably  from  the  polarization  energy 


increase . 
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3.  The  (R  H)+  clusters  have  the  structure  of  a  (RH)+  molecule  and  n-1 

n 

neutral  R  atoms.  The  neutral  rare  gas  atoms  are  located  at  the  H  side  of  the 

(ArH)+  molecule  in  the  (ArH)+Ar  .  clusters  and  at  the  Xe  side  of  the  (XeH)  + 

n-  i 

molecule  in  the  (XeH)+Xe  .  clusters. 

n- 1 

4.  In  order  to  consider  the  (R  H„)+  clusters,  we  have  modified  the  DIIS 

n  l 

method  by  including  into  the  DIIS  wave  function  some  DIM  components  which 

describe  ^2  and  states  of  the  R+H  fragments.  We  find  that  the  Xe+  ion 

forms  with  H2  molecule  a  stable  (Dg  -  0.46  eV)  collinear  cluster  (XeHH)+  with 

the  strong  contribution  of  valence  forces.  When  more  Xe  atoms  are  involved, 

the  cluster  is  formed  by  Xe+  (n  >  1)  ion  and  H„  molecule  which  are  bound  one 

n  z 

to  another  by  polarization  forces  only.  In  the  case  of  Ar  atoms,  we  do  not 
find  any  stable  ionic  (ArH2)+  clusters.  The  quasistable  clusters  are  formed 
by  slightly-charged  Ar  atoms  and  ionic  molecule. 
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Appendia ; _ Diatomic  Potentials  and  Exchange  Terms 

The  matrix  elements  within  the  DIIS  method  are  expressed  in  terms  of  the 

potentials  and  exchange  terms  of  diatomic  fragments  (see  Sec.  II).  In  the  R*, 

(R  H)+  and  (R  H„)+  ionic  clusters  considered  in  this  work,  there  are  the 
n  n  z 

+  +  + 

following  neutral  and  ionic  diatomic  fragments:  H2>  RH,  R2>  H2  and  (RH)  . 

The  potentials  of  the  neutral  fragments,  W^ ,  are  used  directly  in  the  matrix 
elements  expressions.  The  adiabatic  potentials  of  the  ionic  fragments. 
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W^j  ,  are  needed  for  the  calculations  of  the  diabatic  potentials  (see  Eqs . 
(3),  (4)  and  (29)).  When  the  heteronuclear  fragment  (RH)  +  is  considered,  one 
needs  also  the  exchange  term  V  which  has  to  be  defined  independently  from  the 
potentials.  The  potentials  are  given  in  eV  and  the  distances  in  A. 


A.  Neutral  Fragments  and  RH 


We  need  the  following  potentials  of  the  neutral  fragments:  1SR£,  12^H, 
13  + 

Zjyj  and  (the  last  one  is  for  (R^i^)  only).  The  potentials  are  taken 


as  a  sum  of  a  repulsive  exponential  term  and  an  attractive  dispersion  term. 


WArAr(ls)  "  6940  exPl-3-6Rl  -  6°/R6  ,  (A.  1) 

68 

which  provides  the  empirical  values  Rg  -  3.76  A  and  Dg  -  0.012  eV  and  fits 

69 

the  afe  initio  potential  at  small  R,  and 


WXeX/S)  -  24000  exp{  -  3 . 17R]  -  320/R6  ,  (A.  2) 

which  provides  the  empirical7*3  values  R^  -  4.47  A  and  Dg  -  0.023  eV  and  fits 

the  afe  initio  potential71  at  small  R. 

The  ground-state  potential  is  described  by  a  modified  (for  large  R) 
58 

Morse  function, 


“hA*  -  4.746  G(G-2)g  (A. 3) 
G  -  exp( -a(R-Re) ]  (A. 4) 
g  -  exp[-b(R-Re)3] 


(A.  5) 
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In  Eqs .  (A. 4)  and  (A. 5),  Rg  -  0.7417  A,  and  the  parameters  a  -  1.9446  A*^  and 

•3  72 

b  -  0.1215  A  are  found  by  fitting  the  exact  ^  potential.  The  accuracy  of 

the  fit  is  generally  within  0.05  eV  for  R  >  0.53  A.  The  repulsive  2 

potential  of  is  described  by  an  anti -Morse  function, 


Wjjjj CZ)  -  1.963  G(G-2)g  , 


(A. 6) 


where  the  G-function  parameter  is  a  -  1.805  A and  the  g-function  parameter 

-  3  72 

is  b  -  0.1215  A  The  potential  (A. 6)  fits  the  exact  potential. 

The  RH  potential  is  described  by  a  Lennard- Jones  potential  with  the 

7  3 

empirical  parameters  for  ArH  of  Rg  -  3.62  A  and  Dg  -  4.15  meV  and  for 


XeH  of  R  -  3.93  A  and  D  -  6.85  meV: 
e  e 
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WArH  (ls)  “  O-OlfiSU^i2)12  *  <**£*> *1 


WXeH(1S)  -  0.0274[  <^) 12  -  (^)6] 


(A.  7) 


(A.  8) 


B.  Ionic  Homonuclear  Fragments  R*  and 

We  need  the  following  adiabatic  potentials  of  the  ionic  homonuclear 

2  2  2  2  +  +  2 
fragments :  2  ,  E  ,  n  ,  n  for  R„  and  a  ,  a  for  H„ .  The  attractive  E 

6  u’  g’  g’  u  2  g’  u  2  u 

2  + 

and  H  potentials  of  R«  are  described  as  a  sum  of  Morse  function  and  the 
S 

polarization  term  P.  For  argon  these  are  expressed  as 


W/4  a  x+(2S  )  -  1.184  G(G- 2)  +  PA 
(ArAr)  u'  '  Ar 

W/4  *  %+(2n  )  -  0.025  G(G-2)  +  PA 
(ArAr)  g  Ar 


R  -  2.48  ,  a  -  2.0  (A. 9) 

e 

R  -  3.2  ,  a  -  2.3  (A. 10) 

e 
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PAr  -  -11.7(1  +  (^)12]'1/3/R4  .  (A. 11) 

where  Rg  and  a  in  (A. 9)  and  (A. 10)  are  the  parameters  of  the  G-function  (A. 4), 
and  the  number  3.8  in  (A.  11)  is  twice  the  radius  of  the  Ar  atom.  We  assume 
the  Ar+  radius  to  be  equal  to  the  neutral  Ar  radius  rAr  -  1.9  A.  For  xenon  we 
have 

W/v  v  ,+(2Z  )  -  1.0  G(G-2)  +  Pv  ,  R  -  3.22  ,  a  -  155  (A. 12) 

(X.exe)  u  xe  e 

W.v  v  .+(2n  )  -  0.01  G(G-2)  +  Pv  .  R  -  3.8  ,  a  -  3.0  (A. 13) 

(.AeXe/  g  Ae  e 

PXe  -  -2.91(1  +  (^)12]*1/3/R4  •  (A. 14) 

The  Xe  and  Xe+  radii  are  both  taken  as  equal  to  r^  -  2.2  A,  so  that  twice  the 
double  radius  in  (A. 14)  is  4.4  A. 

The  repulsive  potentials  of  R^  are  described  by  a  sum  of  a  polarization 

term  P  and  anti -Morse  function  for  the  Z  state  or  exponential  function  for 

g 

the  II  state: 
u 


W(ArAr)+<2V  "  0  9  G(G+2)  +  PAr  ’  Re  -  2.48  ,  a  -  1.72  (A.  15) 

B<ArAr>+<\>  *  1290  «*W-2.Sal  *  PAr  (A.  16) 

W,v  _  ,+(2Z  )  -  0.78  G(G+2)  +  P„  ,  R  -  3.22  ,  a  -  1.45  (A. 17) 

(AM6;  g  Ae  e 

'  6-9l°  «Pt-2  75Rl  *  PXe  .  (A.  18) 
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The  ground-state  potentials  (A. 9)  and  (A. 12)  provide  the  empirical 
dissociation  energies^  of  Dg  -  1.24  eV  for  Arj  and  Dg  -  1.08  eV  for  Xe^.  The 
excited-state  potentials  (A. 10),  (A. 13)  and  (A.  15)  -  (A.  18)  fit  the  ab  initio 
potentials  of  Ref.  75. 

The  exact  ground- state  potential  a  of  the  molecular  ion^  has  been 
fit  by  a  Morse  function  in  Ref.  77  (there  is  a  misprint  in  Ref.  77  on  p.  17: 
a  -  0.6678  should  be  a  -  0.678).  In  order  to  improve  the  fitting  at  small  R, 
we  have  added  a  repulsive  exponential  component  to  the  Morse  function.  At 
large  R,  the  fitting  is  improved  by  the  g-function  (A. 5), 

W  (a  )  -  810.7  exp[ -11.34R]  +  2.795  G(G-2)g  (A. 19) 

(HH)  g 

with  the  following  parameters  for  the  G-  and  g-functions:  R  -  1.054  A,  a  - 
1.281  A,  b  -  0.0067  A"3.  The  accuracy  of  the  fit  is  within  0.015  eV  for  R  > 
0.74  A.  The  equilibrium  distance  of  the  potential  (A.  19)  is  1.06  A,  slightly 
larger  than  the  parameter  R@  for  the  G-  and  g-functions.  The  dissociation 
energy  of  the  potential  (A.  19)  is  Dg  -  2.79  eV.  The  excited  a u  potential ^  is 
fit  by  an  anti -Morse  potential, 

W  (a  )  -  3.014  G(G+2)g  ,  R  -  1.054  ,  a  -  1.379  ,  b  -  0.0067  . 

(HH)  U  * 

(A. 20) 


(A. 21) 
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where  1.7  is  the  H  radius . 


C.  Ionic  Heteronuclear  Fragment  + 

In  order  to  find  from  Eq.  (4)  or  (29)  the  (RH)+  diabatic  Z  potentials  U, 
one  needs  three  adiabatic  potentials,  namely  ^Z  and  3Z  asymptotic  to  R+  +  H 
and  Z  asymptotic  to  R  +  H  .  The  adiabatic  potentials  are  described  as  a  sum 
of  an  attractive  (Morse)  or  repulsive  (anti-Morse  or  exponential)  potential 
and  the  polarization  term  P .  The  adiabatic  Z  potentials  of  (ArH)+  are 
expressed  as  follows: 


W(ArH)+<xls)  “  3‘157  G(G‘2)  +  pAr  -  R#  -  1.266  ,  a 
W(ArH)+(Bl2)  “  2'75  G(G+2)  +  PH  +  2‘16  •  Re  -  1-266 

W(ArH)+<3S)  "  600  exp[-1.8R]  +  P^  +  2.16 
PAr  -  11.7(1  +  (J^2)12]'1/3/R4 
PH  -  4.8(1  +  (^)12]'1/3/R4  . 


1.94  (A. 22) 

a  -  1.2  (A. 23) 

(A.  24) 

(A. 25) 

(A. 26) 


In  (A. 22)  and  (A. 23),  Rg  and  a  are  the  parameters  for  the  G-function  (A. 4), 
and  in  (A. 23)  and  (A. 24)  the  number  2.16  presents  the  difference  between  the 
ionization  potentials  of  Ar  and  H.  In  the  polarization  terms  (A. 25)  and 
(A. 26),  1.9  is  the  Ar  radius  and  3.6  is  the  sum  of  the  Ar  and  H  radii.  The 
ground-state  (ArH)+  equilibrium  distance  Rg  -  1.266  k  in  the  potential  (A.  22) 
is  equal  to  the  experimental  value .  ’  This  potential  provides  the  ajjt 


34 


80 

initio  dissociation  energy  Dg  -  4.055  eV.  The  &&  initio  ground- state 

potential  of  Ref.  80  is  used  to  find  the  parameter  a  of  (A. 22).  The  excited- 

state  potentials  (A. 23)  and  (A. 24)  are  found  by  fitting  the  all  initio 

potentials  of  Ref.  81.  The  (ArH)+  exchange  term  V  is  found  from  Eq.  (4)  by 

80 

fitting  the  all  initio  dipole  moment,  which  is  close  to  the  experimental 
,  81 

value  at  the  equilibrium  point, 

V  -  -95.5  R2  exp[-2.42R]  .  (A.27) 

The  adiabatic  Z  potentials  of  (XeH)+  are  as  follows: 

W(XeH)+<xl2)  -  4.03  G(G-2)  +  ?H  ,  Rg  -  1.61  ,  a  -  1.56  (A. 28) 

W(XeH)+<Bl^  “  73  ®xp[ -1.3R)  +  PXe  +  1.47  (A. 29) 

W(XeH)+<32)  “  1.3  G(G+2)  +  PH  ,  Re  -  1.61  ,  a  -  1.56  (A. 30) 

PR  -  -4. 8[1  +  (^)12]'1/3/R4  (A. 31) 

PXe  -  -2.91[1  +  (2^2)12]'1/3/R4  .  (A. 32) 

In  (A. 29),  1.47  presents  the  difference  between  the  ionization  potentials  of  H 
and  Xe.  In  the  polarization  terms  (A. 31)  and  (A. 32),  3.9  is  the  sum  of  the  Xe 
and  H  radii  and  2.2  is  the  Xe  radius.  The  ground-state  (XeH)+  equilibrium 

distance  R@  -  1.61  A.  in  the  potential  (A. 28)  is  equal  to  the  experimental 

82  85 

value.  The  potential  (A. 28)  provides  the  sfe  initio  dissociation  energy  Dg 

-  4.05  eV,  a  value  close  to  the  experimental  dissociation  energy  of  Dg  »  3.8 
84 

eV .  The  parameter  a  of  (A.  28)  is  found  by  fitting  the  gli  initio  potential 
of  Ref.  83.  The  excited-state  potentials  (A. 29)  and  (A. 30)  are  found  by 


fitting  the  ab  initio  potentials  of  Ref.  85.  The  (XeH)+  exchange  term  V  is 

83 

found  from  Eq.  (4)  by  fitting  the  at£  initio  dipole  moment, 

V  -  -57.7  R2  exp [-2.0  R]  .  (A. 33) 

The  II  potentials  of  (ArH)+  are  obtained  by  fitting  the  ak  initio 
potentials  of  Ref.  81, 

W(ArH)+(AllI)  “  423  e*P(-3-9Rl  +  +  2.16  (A. 34) 

W(ArH)+(3lI)  “  90  ®xp[-1.2R]  +  ?H  +  2-16  ■  (A-35> 

where  P is  given  by  (A. 26).  The  II  potentials  of  (XeH)+  are  obtained  by 
fitting  the  afe.  initio  potentials  of  Ref.  85, 

W/v  ux+A  -  370  exp[  -  2 . 74R)  +  (A.  36) 

W(XeH)+(3lI)  "  430  exp[-2.7^RJ  +  ?H  ,  (A. 37) 

where  P„  is  given  by  (A. 31). 
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Table  1.  Structure  of  the  ionic  clusters  Ar+  and  ionic  molecule  Ar*.  The 

n  2 

distances  (R)  are  in  A,  and  the  energies  (E,D)  in  eV. 


N 

Ar+ 

n 

Structure  ^ 

r2) 

r3> 

q4> 

E 

Dissociation 

Products 

D 

I 

Ar2 

2.48 

-1.240 

Ar*  +  Ar 

1.240 

II 

at; 

linear 

2.59 

-1.443 

Ar*  +  Ar 

0.203 

III 

Ar4 

(Ar*)  (Ar) 

3  z  xz 

2.59 

3.68 

0.002 

-1.490 

Ar*  +  Ar 

0.047 

IV 

A< 

regular  pyramid 

2.836 

-1.480 

Ar*  +  2Ar 

0.240 

V 

A< 

linear 

2.59 

3.3 

0.01 

-1.474 

Ar*  +  Ar 

0.031 

VI 

Ar5 

<Ar3>2<Ar2>xy 

2.59 

3.68 

0 

-1.540 

(Ar4*IH  + 

Ar 

0.050 

VII 

Ar5 

<Ar3>z(Ar2>x2 

2.59 

3.61 

0.002 

-1.540 

(Ar4>III  + 

Ar 

0.050 

VIII 

Ar* 

<Ar3>z(A'2>xz 

2.59 

3.68 

0.003 

-1.537 

<Art>III  + 

Ar 

0.047 

IX 

Ar! 

<Ar3VAt2VAr)x 

2.59 

3.61 

0.001 

-1.601 

(Ar*)VI  +  Ar 

0.061 

^The  subscripts  z,  xz,  etc.  denote  the  location  along  an  axis  or  in  a  plane  (see 
Fig.  l.II-V). 

2) 

Distance  between  adjacent  charged  atoms. 

3) 

Distance  from  a  neutral  or  weakly -charged  atom  to  the  nearest  charged  atom. 

4) 

Common  charge  of  all  weakly -charged  atoms. 


41 


Table  2.  Transition  energy  T  in  eV  and  transition  moment  n  in  Debyes 
of  the  allowed  transitions  in  Ar^  and  Ar+.  q  is  the  charge 
distribution  in  the  excited  states. 


nL) 

Ar+ 

n 

T 

D 

1 

q 

2 

3 

4 

II 

Ar3 

2.26 

8.8 

0.5 

0 

0.5 

III 

Ar*Ar 

2.24 

6.8 

0.30 

0.01 

0.34 

0.35 

2.19 

2.3 

0.10 

0.04 

0.24 

0.62 

1.79 

0.7 

0.27 

0.37 

0.15 

0.21 

V 

Art  linear 

4 

1.97 

8.2 

0.24 

0.01 

0.26 

0.49 

IV 

Ar*  pyramid 

1.46 

2.1 

0.25 

0.25 

0.25 

0.25 

1.46 

2.1 

0.30 

0.30 

0.20 

0.20 

1.46 

2.1 

0.20 

0.20 

0.30 

0.30 

^For  N,  see  Table  1. 
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Table  3.  Structure  of  the  ionic  clusters  Ar^(ArH)+  and  ionic  molecule  (ArH)+. 

The  distances  (R)  are  in  A,  the  energies  (E)  in  eV  and  the  dipole 
moment  (/i)  in  Debyes. 


Ark(ArH)+ 

R  A  u 

Ar-H 

2) 

ra  * 

Ar-Ar 

2) 

M 

E 

Ar 

Detachmei 

(ArH)+ 

0.569 

2.78 

-4.055 

Ar(ArH)+ 

0.579 

2.81 

3.70 

5.98 

-4.218 

0.163 

Ar2(ArH)+ 

0.591 

2.79 

3.81 

6.32 

-4.395 

0.177 

Ar3(ArH)+ 

0.601 

2.79 

3.72 

6.23 

-4.586 

0.191 

Ar4(ArH)+ 

0.609 

2.84 

3.57 

4.94 

-4.760 

0.174 

Arg(ArH)+ 

0.614 

2.84 

3.17 

3.87 

3.38 

3.33 

-4.947 

0.094 

^ Charge  of  the  H  atom  in  (ArH)+ 

2) 

R.  ,,  and  R.  .  are  the  distances  from  the  neutral  atom  to  H  and  Ar  of 
Ar-H  Ar-Ar 

(ArH)+,  respectively.  For  Ar,(ArH)+  there  are  two  different  Ar-H  and  Ar-Ar 

o 


distances . 
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Table  4.  Structure  of  the  ionic  clusters  Ar  Ht  and  ionic  molecule  Ht.  The 

n  l 

distances  (R)  are  in  A,  the  energies  in  eV  and  the  dipole  moments 
(n)  in  Debyes. 


ArnH2 

qH+ 

H2 

^-H 

X) 

Ar 

M 

E 

Ar 

Detachment 

T 

Transitions 
_ a _  T 

“2 

1.0 

1.06 

0 

-2.790 

11.8 

2.5 

ArH^ 

0.98 

1.05 

2.935 

12.3 

-2.983 

0.193 

5.65 

1.9 

Ar2H2 

0.964 

1.05 

2.94 

9.02 

-3.183 

0.200 

5.83 

2.0 

Ar3H2 

0.95 

1.04 

2.95 

5.15 

-3.377 

0.194 

5.58 

0.8 

5.89 

2.3 

Ar4H2 

0.936 

1.03 

2.955 

1.65 

-3.568 

0.191 

5.71 

0.9 

5.94 

2.5 

ArX 

0.926 

1.03 

2.995 

0 

-3.754 

0.186 

5.44 

1.7 

5.92 

2.2 

13  + 

‘'Distance  from  the  Ar  atoms  to  the  center  of  The  adjacent  Ar  atoms  are 

separated  by  3.7  A  (n  -  2-4)  and  3.52  A  in  symmetrical  Ar^H^- 
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Table  5.  Structure  of  the  ionic  clusters  Xe+  and  ionic  molecule  Xet.  (See 

n  2. 

Table  1  for  the  footnotes . ) 


N 

Xe+ 

n 

structure  ^ 

r2> 

r3> 

<4) 

E 

Dissociation 

Products 

D 

I 

Xe2 

3.22 

-1.080 

Xe'1 

'  +  Xe 

1.08 

II 

Xe3 

linear 

3.383 

-1.277 

Xe2 

+  Xe 

0.197 

III 

XeI 

regular  pyramid  3 . 646 

-1.480 

Xe* 

+  2Xe 

0.40 

IV 

X< 

(Xe3>z<Xe)xz 

3.384 

4.19 

0.014 

-1.356 

Xe3 

+  Xe 

0.079 

V 

X< 

linear 

3.315 

3.65 

0.184 

-1.346 

Xe* 

+  Xe 

0.069 

VI 

Xe5 

K>III(Xe)xz 

3.646 

4.40 

0.001 

-1.568 

X< 

+  Xe 

0.088 

VII 

Xe5 

<XetVXVxz 

3.383 

4.24 

0.005 

-1.439 

Xe*Xe 

+  Xe 

0.083 

VIII 

Xe6 

^Xe4^IIIXe2 

3.644 

4.40 

0 

-1.637 

Xe*Xe 

+  Xe 

0.069 

IX 

Xel 

(Xe3>z^Xe3^xz 

3.375 

4.20 

0.021 

-1.547 

Xe^Xe2  +  Xe 

0.108 
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Table  6.  Transition  energies  T  (eV)  and  moments  m(D)  of  the  Xe^Xe^  allowed 
transitions. 


N  X> 

XetXe, 

3  k 

T 

T 

Transitions 

M  T 

II 

Xe* 

1.87 

11 

IV 

Xe*Xe 

1.86 

8.3 

1.66 

3.6 

1.17 

1.7 

VII 

Xe^Xe2 

1.86 

11 

1.15 

1.7 

IX 

Xe^Xe  3 

1.88 

7.9 

1.82 

4.5 

1.23 

1.6 

^For  N,  see  Table  6. 
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Table  7. 

Structure  of 

the  ionic 

clusters  Xe^CXeH) 

+  and 

ionic  molecule  (XeH) 

Xek(XeH)+ 

D  +  2) 

nCe-Xe 

M 

(V) 

E 

Xe 

Detachment 

(XeH)+ 

0.282 

1.54 

-4.050 

Xe(XeH)+ 

0.280 

6.07 

4.46 

11.5 

-4.096 

0.046 

Xe2(XeH)+ 

0.279 

5.95 

4.50 

12.7 

-4.156 

0.060 

Xe3(XeH)+ 

0.278 

5.96 

4.56 

13.3 

-4.235 

0.079 

Xe4(XeH)+ 

0.277 

5.91 

4.62 

12.4 

-4.297 

0.062 

Xe5(XeH)+ 

0.275 

5.73 

4.85 

10.7 

-4.372 

0.038 

^ Charge  of  the  H  atom  In  (XeH)+. 

2) 

Ae  H  an<1  RXe  Xe+  are  t*ie  distances  from  the  neutral  atom  to  H  and  Xe  of 
(XeH)+,  respectively. 
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Figure 

Fig.  1 

Fig.  2 

Fig.  3 

Fig.  4 


Captions 

Ar ^  and  Ar+  geometry.  The  Roman  numerals  indicate  the  clusters  (the 
same  as  in  Table  1).  The  numbers  without  a  sign  are  distances  in  k, 
and  the  numbers  with  the  +  sign  are  the  atomic  charges. 

Geometry  of  the  ionic  clusters  Ar(ArH)+,  Ar£(ArH)+  and  ArH^.  For 
notations,  see  Fig.  1. 

Geometry  of  Xe+  and  Xe*.  For  notations,  see  Fig.  1.  The  Roman 
numerals  indicate  the  clusters  (see  Table  5). 

Geometry  of  the  ionic  clusters  Xe(XeH)+,  Xe2(XeH)+,  (Xel^)*  and 
Xe^H^.  For  notations,  see  Figs.  1  and  3. 
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